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We present a new semi-empirical model for calculating electron transport in atomic-scale devices. 
The model is an extension of the Extended Hiickel method with a self-consistent Hartree potential. 
This potential models the effect of an external bias and corresponding charge re- arrangements 
in the device. It is also possible to include the effect of external gate potentials and continuum 
dielectric regions in the device. The model is used to study the electron transport through an organic 
molecule between gold surfaces, and it is demonstrated that the results are in closer agreement with 
experiments than ab initio approaches provide. In another example, we study the transition from 
tunneling to thermionic emission in a transistor structure based on graphene nanoribbons. 
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INTRODUCTION 



are their lower computational cost. 



As the minimum feature sizes of electronic devices are 
approaching the atomic scale, it becomes increasingly 
important to include the effects of single atoms in de- 
vice simulations. In recent years, there have been sev- 
eral developments of atomic-scale electron transport sim- 
ulation models based on the Non-Equilibrium Green's 
Function (NEGF) formalism^. The approaches can 
roughly be divided into two catagories: ab initio ap- 
proaches, where the electronic structure of the system 
is calculated from first principles, typically with Den- 
sity Functional Theory (DFT)^ ^, and semi-empirical ap- 
proaches, where the electronic structure is calculated us- 
ing a model with adjustable parameters fitted to experi- 
ments or first-principles calculations. Examples of semi- 
empirical transport models are methods based on Slater- 
Koster tight-binding parameters^ ^ and Extended Hiickel 
parameter^^HEl. 

The ab initio models have the advantage of predictive 
power, and can often give reasonable results for systems 
where there is no prior experimental data. However, the 
use of the Kohn-Sham one-particle states as quasiparti- 
cles is questionable, and it is well known that for many 
systems the energies of the unoccupied levels are rather 
poorly described within DFT. Furthermore, solving the 
Kohn-Sham equations can be computationally demand- 
ing, and solving for device structures with thousands of 
atoms is only feasible on large parallel computers. 

The semi-empirical models have less predictive power, 
but when used within their application domain they can 
give very accurate results. The models may also be fit- 
ted to experimental data, and can thus in some cases give 
more accurate results than DFT-based methods. How- 
ever, the main advantage of the semi-empirical methods 



In this paper we will present the formalism behind a 
new semi-empirical transport model based on the Ex- 
tended Hiickel (EH) method. The model can be viewed as 
an extension to the work by Zahid et al}'^ ^ with the main 
difference being the treatment of the electrostatic interac- 
tions. Zahid et al. only describe part of the electrostatic 
interactions in the device; most importantly, they use 
the Fermi level of the electrodes as a fitting parameter 
and do not account for the charge transfer from the elec- 
trodes to the device. In the current work, the Fermi level 
of the electrodes is determined self-consistently by using 
the methodology introduced by Brandbyge et al.^^ln this 
way, we include the charge transfer from the electrodes 
to the device region and describe all electrostatic interac- 
tions self-consistently. This is accomplished by defining 
a real-space electron density and numerically solving for 
the Hartree potential on a real-space grid. Through a 
multi-grid Poisson solver, we include the self-consistent 
field from an applied bias, and allow for including con- 
tinuum dielectric regions and electrostatic gates within 
the scattering region. 



The organization of the paper is the following: In sec- 
tion [IT] we introduce the self- consistent Extended Hiickel 
(EH-SCF) model, and in section 



HI 



we present the for- 
malism for modelling nano-scale devices. In section IV 
we apply the model to a molecular device, and in sec- 
tion [V| we consider a graphene nano-transistor where an 
electrostatic gate is controlling the electron transport in 



the device. Finally, in section VI, we conclude the paper. 
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II. THE SELF-CONSISTENT 
EXTENDED-HUCKEL METHOD 



following approximation for the spatial distribution of the 
induced electron density: 



In this section we describe the EH-SCF framework. In 
Extended Hiickel theory, the electronic structure of the 
system is expanded in a basis set of local atomic orbitals 
(LCAOs) 



(1) 



where Yj^ is a spherical harmonic and Rni is a superpo- 
sition of Slater orbitals 



Rni{r) 



(2n)! 



(2) 

The LCAOs are described by the adjustable parameters 
^1 7 ^2 , C'l , and C2 , and these parameters must be defined 
for the valence orbitals of each element. 

The central object in EH theory is the overlap matrix, 

C _ / ^ij ~ 

' ~\Iv - - Rj) dr if R, ^ R, 

(3) 

where Hs a composite index for nlm and R^ is the posi- 
tion of the center of orbital i. 

From the overlap matrix, the one-electron Hamiltonian 
is defined by 



Sn{r) = Sm^ 



(6) 



of each Gaussian is 



where the weight (5m^ = 
the excess charge of atom ja as obtained from the Mul- 
liken population and the ion valence charge Z^. Sub- 
sequently, the Hartree potential is calculated from the 
Poisson equation 



V- [e(r) VSVH{r)] = Sn{r), 



(7) 



(4) 

where Ei is an orbital energy, and Pi is an adjustable pa- 
rameter (often chosen to be 1.75). SVH(R^i) is the Hartree 
potential corresponding to the induced electron density 
on the atoms, i.e. the change in electron density com- 
pared to a superposition of neutral atomic-like densities. 
This term must be determined self-consistently, and is 
not included in standard EH modelJ^^. In the following 
section we describe how this term is calculated. 



which is solved with the appropriate boundary conditions 
on the leads and gate electrodes imposed by the applied 
voltages. Here, e(r) is the spatially dependent dielectric 
constant of the device constituents, and allows for the 
inclusion of dielectric screening regions. 

To see the significance of the width of the Gaussian 
orbital, let us calculate the electrostatic potential from 
a single Gaussian electron density at position R^. The 
result is 

Erff^/c^lr-R.!) 
5VH{r) = e(m, - Z,) (8) 

and from this equation we see that the on-site value of 
the Hartree potential is 5VH{^fi) = (^/i — Z^)^^^ where 
the parameter 
\ii= j 

S^3 if^y^J ^^ = 2eJ^ (9) 



is the on-site Hartree shift. The param eter 7^ is a well- 
known quantity in CNDO theorj^i^^, and values of 7^ 
are listed for many elements in the periodic table. Thus, 
we fix the value of 7^ using CNDO theory, and then use 
Eq. (|9| to calculate the value of for each element. 



HI. EH-SCF METHOD FOR A NANO-SCALE 
DEVICE 



A. Solving the Poisson Equation to Obtain the 
Hartree Potential 

To calculate the induced Hartree potential we need to 
determine the spatial distribution of the electron den- 
sity. To this end, we introduce the Mulliken population 
of atom number ji 

where Dij is the density matrix. The total number of 
electrons can now be written as a sum of atomic contri- 
butions, N = m^. 

We will represent the spatial distribution of each 
atomic contribution by a Gaussian function, and use the 



Fig. [T] illustrates the setup of a molecular device sys- 
tem. The system consists of three regions: the central 
region, and the left and right electrode regions. The cen- 
tral region includes the active parts of the device and 
sufficient parts of the contacts, such that the properties 
of the electrode regions can be described as bulk materi- 
als. For metallic contacts, this will typically be achieved 
by extending the central region 5-10 A into the contacts. 

The calculation of the electron transport properties of 
the system is divided into two parts. The first part is 
a self-consistent calculation for the electrodes, with peri- 
odic boundary conditions in the transport direction. In 
the directions perpendicular to the transport direction, 
we apply the same boundary conditions for the two elec- 
trodes and the central region, and these boundary con- 
ditions are described below. 
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FIG. 1. (Color online) Geometry of a nano-device consist- 
ing of a dithiol-triethynylene-phenylene molecule attached to 
two (3x3) (111) gold electrodes. The left and right electrode 
regions are illustrated with wire boxes, and the properties 
of these regions are obtained from a calculation with periodic 
boundary conditions in all directions. The region between the 
two electrodes is the central device region, which is described 
with open boundary conditions in the transport direction, and 
periodic boundary conditions in the directions perpendicular 
to the transport direction. 



In the second part of the calculation, the electrodes 
define the boundary conditions for a self-consistent open 
boundary calculation of the properties of the central re- 
gion. The main steps in the open boundary calculation is 
the determination of the density matrix, the evaluation 
of the real-space density, and, finally, the calculation of 
the Hartree potential. These steps will be described in 
more detail in the following section. 



A. Calculating the Self-Consistent Density Matrix 
of the Central Region 

In this section we will describe the calculation of the 
density matrix of the central region. We assume that the 
self-consistent properties of the left and right electrodes 
have already been obtained, and thus we also know their 
respective Fermi levels, and e^. We allow for an exter- 
nal bias Vh to be applied between the two electrodes, and 
define the left and right chemical potentials /iL = 
and ijlr = e^. The applied bias thus shifts all energies 
in the left electrode, and a positive bias gives rise to an 
electrical current from left to right. 

The density matrix for this non-equilibrium system, 
with two different chemical potentials, is found by filling 
up the left and right originating states according to their 
respective chemical potentials^ 



by the self energies and S^, which arise due to the 
coupling of the central region with the semi-infinite left 
and right electrodes, respectively. 

Further details of the NEGF formalism can be found 
in Refs. 14, Here we just note that to improve the nu- 
merical efficiency, the integral in Eq. (10) is divided into 



an equilibrium and non-equilibrium part. The equilib- 
rium part is calculated on a complex contour far from the 
real-axis poles of the Green's function, and only the non- 
equilibrium part is performed along the real axis. The 
equilibrium and non-equilibrium parts are then joined 
using the double-contour technique introduced by Brand- 
byge et al^ 

From the density matrix we may now evaluate the real- 
space density in the central region using Eq. (|6|. It is 
important to note that near the left and right faces of 
the central region there will be contributions from the 
electrode regions, and this "spill in" must be properly 
accounted for. 

Once the real-space density is known, the Hartree po- 
tential is calculated by solving the Poisson equation in 
Eq. using a real-space multi-grid method. On the left 
and right faces of the central region the Hartree potential 
is fixed by the electrode Hartree potentials, appropriately 
shifted according to the applied bias. In the directions 
perpendicular to the transport directions, we apply the 
appropriate boundary conditions, fixed or periodic, as 
demanded by e.g. the presence of gate electrodes. 

The so-obtained Hartree potential defines a new 
Hamiltonian, via Eq. [4j and the steps in section HI A 
must be repeated until a self-consistent solution is ob- 
tained. 



B. Transmission and Current 

Once the self-consistent one-electron Hamiltonian has 
been obtained, we can finally evaluate the transmission 
coefficients^ 

T{e) = Tr[tde)G\e)tR{e)G{e)] (12) 
and the current 



D 



2e 

[p''{e)nF{e - ^l) ^ p'^{e)nF{e - m)] de, j _J^(^)[^F{e - pl) - np^e - pR)]de. (13) 



(10) 

where {p^) is the contribution to the spectral den- 
sity of states from scattering states originating in the left 
(right) reservoir. 

The calculation of the spectral densities is performed 
using NEGF theory, and we write the partial spectral 
densities as 

pL,R^e) = ^G{e)t^-^{e)GHe), (11) 

where G is the retarded Green's function of the central 
region, and the broadening function F = — S'''] is given 



In the following sections, we apply this formalism to 
the calculation of the electrical properties of a molecule 
between gold electrodes, as well as a graphene nano- 
transistor. 



IV. TOUR WIRE BETWEEN GOLD 
ELECTRODES 

In this section we will investigate the electrical prop- 
erties of a phenylene ethynylene oligomer, also popularly 
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called a Tour wire. We will compare the electrical proper- 
ties of the molecule when it is symmetrically and asym- 
metrically coupled with two Au(lll) surfaces. In the 
symmetric system, as illustrated in Fig. [l] the molecule 
is connected with both gold surfaces through thiol bonds, 
whereas the asymmetric system only has a thiol bond to 
one of them. 

The system has previously been investigated exper- 
imentally b y K ushmerick et al^^ and theoretically by 
Taylor et a/. , and it has been found that the asymmet- 
rically coupled system shows strongly asymmetric I-V 
characteristics^^ . 

The calculations by Taylor et al. were based on DFT- 
LDA, and the asymmetric behaviour could be related to 
the voltage drop in the system. This system is therefore 
an excellent testing ground for our semi-empirical model, 
since a correct description of the electrical properties re- 
quires not only a good model for the zero-bias electronic 
structure, but also a good description of the bias-induced 
effects. 



A. Transmission Spectrum of the Symmetric Tour 
Wire Junction 

To setup the symmetric system we first relaxed the 
isolated Tour wire using DFT-LDA^. During the re- 
laxation, passivating hydrogen atoms were kept on the 
sulfur atoms. Afterwards, these hydrogen atoms were re- 
moved and the two sulfur atoms placed at the FCC sites 
of two Au(lll)-(3x3) surfaces. The height of the S atom 
above the surface was 1.9 A (corresponding to an Au-S 
distance of 2.53 A). 

We nex t set up the EH model with Hoffmann 
parameters"*^^'^^ and perform a self-consistent calculation 
to obtain transmission spectra for different k-point sam- 
pling grids. The results are shown in the upper plot of 
Fig. [2] In each case, the same k-point grid was used for 
both the self-consistent and transmission calculation, and 
we see from the figure that using (1x1) k-point is insuffi- 
cient while (2x2) and (4x4) k-points give almost identical 
results. Thus, we will use a (2x2) k-point sampling grid 
for the remainder of this study. 

In the lower plot of Fig. [2] we compare the transmis- 
sion spectra calculated with DFT-LDA, EH-SCF, and 
EH without the Hartree term of Eq. Q. For the DFT- 
LDA model we use similar parameters as Taylor et al^^ 
except for the k-point sampling which is (2x2) in the cur- 
rent study. The calculations in Ref. "21" were performed 
with a (1x1) k-point sampling, which is insufficient^^, and 
thus the DFT-LDA results in this study will differ from 
those by Taylor et al. 

For the EH calculation we see a peak in the transmis- 
sion spectrum just around the Fermi level of the gold 
electrodes. This peak arises from transmission through 
the LUMO orbital of the Tour wire. 

In the self-consistent EH calculation there will be a 
charge transfer from the gold surface to the LUMO or- 




E-E^^ (eV) 



FIG. 2. (Color online) The upper plot shows the transmis- 
sion spectrum of the symmetric Tour wire device, calculated 
with the EH-SCF model for three different k-point samplings. 
The lower plot shows transmission spectra calculated with 
a (2x2) k-point sampling using different models: EH-SCF 
(solid), EH without the Hartree term (dotted), and DFT- 
LDA (blue dashed). Energies are given relative to the Fermi 
level of the gold electrodes. 



bital, and we see that this gives rise to a shift of the 
orbital by 1 eV, illustrated by the arrow in Fig. [2] 

For the DFT-LDA calculation we see that the LUMO 
peak is shifted further away from the gold Fermi level, 
and the HOMO and LUMO peaks of the transmission 
spectrum are placed almost symmetrically around the 
gold Fermi level. We also note that the transmission co- 
efficient at the Fermi level, corresponding to the zero-bias 
conductance, is almost one order of magnitude higher 
within the DFT-LDA model. We will discuss this fur- 
ther below. 

We also note that Taylor et al. find a LUMO level even 
further away from the gold Fermi level; this is related to 
the insufficient k-point sampling. 



B. I— V Characteristics of the Symmetric and 
Asymmetric Tour Wire Systems 

We will now study both the symmetric and asymmet- 
ric Tour wire system and compare their respective I-V 
characteristics. The geometry of the asymmetric system 
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FIG. 3. (Color online) Geometry of the asymmetric system. 
The Tour wire is attached to the left gold electrode through 
a thiol bond, while the right end of the molecule is hydrogen- 
terminated and there is no chemical bond to the right gold 
electrode. 
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FIG. 4. (Color online) I-V characteristics of the symmetric 
(upper figure) and asymmetric (lower figure) Tour wire device. 
The positive current direction is from left to right. 



is illustrated in Fig. [Sj The geometry is similar to that 
of Fig. [l] except for the right-most sulfur atom which 
has been replaced by a hydrogen atom with a C-H bond 
length of 1.1 A. The distance between the hydrogen atom 
and the right gold surface is 1.5 A. 

We perform self-consistent calculations for both the 
symmetric and asymmetric systems with the EH-SCF 
and DFT-LDA methods, and vary the bias from -1 to 
-hi V in steps of 0.1 V. The results are shown in Fig. |4] 
For the symmetric device we obtain rather similar, sym- 
metric I-V characteristics for both the EH-SCF and 
DFT-LDA methods. The main difference is that the zero- 
bias conductance is significantly higher with DFT-LDA, 
reflecting the higher transmission coefficient at the Fermi 
level, as shown in Fig. [2] 

For the asymmetric device we see that both the DFT- 
LDA and EH-SCF models give rise to rectification - how- 




Au Au Au s cc a: c c oc cc c c oc o: S Au Au Au 
Transport direction 



FIG. 5. (Color online) Voltage drop of the symmetric and 
asymmetric Tour wire systems along a line that goes through 
the two sulfur atoms in the symmetric system, for an ap- 
plied bias of +1 V. The inset shows the voltage drop in the 
asymmetric system subtracted from the voltage drop in the 
symmetric system. Both plots show results calculated with 
the EH-SCF (solid) and DFT-LDA (dashed) models. 



ever, in opposite directions. Taylor et a/.l^ demonstrated 
that the rectification was related to the voltage drop in 
the system, and we therefore in Fig. [5] compare the volt- 
age drops obtained with the two methods. The EH-SCF 
voltage drop is smooth, since the charge density is com- 
posed of a superposition of single rather broad Gaussians 
on each atom. The DFT-LDA model shows atomic-scale 
details, however, as illustrated by the inset, the relative 
change in the voltage drop between the asymmetric and 
symmetric system is quite similar for the EH-SCF and 
DFT-LDA methods. Both methods reveal that in the 
asymmetric system there is an additional voltage drop at 
the contact with the weak bond. This is also one of the 
main results of Taylor et al?'^ . 

The additional voltage drop at the weak contact means 
that the molecular levels of the Tour wire mainly fol- 
low the electrochemical potential of the right electrod^^. 
Since the voltage drop is similar for the DFT-LDA and 
EH-SCF models, the difference in the I-V characteris- 
tics must be related to the different electronic structure 
at zero bias in the two models. Within the DFT-LDA 
model, the transport at the Fermi level is dominated by 
the HOMO. At negative bias, the left electrode has a 
higher electrochemical potential, and electrons from the 
occupied HOMO level can propagate to empty states in 
the left electrode. Thus, for the DFT-LDA model, the 
current is highest for a negative bias at the left electrode. 
For the EH-SCF model, on the other hand, the transport 
at the Fermi level is dominated by the LUMO, and the 
current in this case is highest for a positive bias at the 
left electrode. 

Comparing with the experimental results of Kushmeric 



et 



r 20 



, we find that the EH-SCF rectification direction 
agrees with the experimental rectification direction, while 
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FIG. 6. (Color online) Graphene nano-transistor consisting of 
two metallic zigzag nanoribbons connected by a semiconduct- 
ing armchair ribbon. The nanoribbons are passivated with 
hydrogen, and the width of the ribbons are is 7 A. The device 
is sitting on top of a dielectric and the transport is controlled 
by an electrostatic back-gate. The contour plot illustrates the 
Hartree potential for a gate potential of -1 V. 

the DFT-LDA model predicts rectification in the oppo- 
site direction. We note that the rectification direction 
obtained with our DFT-LDA model is similar to the re- 
sults of Taylor et alW. 

Thus, for this system it seems that the EH-SCF model 
is in better agreement with the experimental results, 
compared to the DFT-LDA model. The example shows 
that the EH-SCF model gives a very good description 
of both the electronic structure and the voltage drop in 
the system. The comparisons between the two methods 
also illustrates how small variations in the positions of 
the HOMO and LUMO levels may change the electrical 
properties of the Tour wire device. 



V. Z-SHAPED GRAPHENE 
NANO-TRANSISTOR 

In this section we will compare the electrical proper- 
ties of a short (34 A) and long (86 A) graphene nano- 
transistor. The system consists of two electrodes consist- 
ing of metallic, zigzag-edge graphene nanoribbons con- 
nected through a semiconducting armchair-edge central 
ribbon. The system is placed 1.4 A above a dielectric 
material with dielectric constant e = 4eo, corresponding 
to Si02. The dielectric is 3 A thick, and below the di- 
electric there is an electrostatic gate. The geometry of 
the short system is illustrated in Fig. [6j A similar system 
was investigated by Yan et al?^ using DFT-LDA. 

For the calculation we use EH parameters from Ref. TS* 
which were derived by fitting to a reference band struc- 
ture of a graphene sheet calculated with DFT-LDA. With 
these parameters, we find a band gap of the central rib- 
bon of 2.2 eV, in agreement with DFT-LDA calculations, 
which illustrates the transferability of the EH parameters 
from 2D graphene to a ID graphene nanoribbon. 



A. Transmission Spectrum 

Fig. [7| shows the transmission spectrum for both the 
long and the short system when there is no applied bias 




E -Ep (eV) 



FIG. 7. Zero-bias transmission spectrum for the short 
(dashed) and long (solid) graphene device. Energies are rela- 
tive to the Fermi level of the electrodes. 

and zero gate potential. The shape of the transmission 
spectrum is directly related to the electronic structure of 
the central semiconducting ribbon. 

The transmission is strongly reduced in the energy re- 
gion from -0.7 to 1.5 eV, corresponding to the band gap 
of the central armchair ribbon. Since there are no energy 
levels in this interval, the electrons must tunnel in order 
to propagate across the junction. For the longer device 
the electrons must tunnel a longer distance, and thus the 
transmission is more strongly reduced. 

Outside the band gap, the transmission is close to 1 
and shows a number of oscillations. Since the central 
ribbon has a finite length, it resembles a molecule with a 
number of discrete energy levels. The levels give rise to 
peaks in the transmission spectrum, and since the longer 
system has more energy levels, the peaks are more closely 
spaced there. 

In the following section we will see how this difference 
in the transmission spectrum gives rise to qualitatively 
different transport mechanisms in the two devices. 



B. Transistor Characteristics 

We now calculate the current for an applied source- 
drain voltage of 0.2 V as a function of the applied gate 
potential. Fig. [8] shows the current for the long and short 
devices, respectively, for gate potentials in the range -1 to 
1 V, for different electrode temperatures. We see that for 
the short device there is only a small effect of the gate 
potential and electron temperature, while for the long 
device the conductance falls off exponentially, reaching a 
minimum in the range to 0.5 V. Moreover, the current 
is strongly temperature-dependent. 

The lack of temperature-dependence for the short de- 
vice shows that the transport is completely dominated 
by electron tunneling. For the long device, on the 
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FIG. 8. (Color online) The tunneling current for a source- 
drain voltage of 0.2 V, as a function of the gate potential 
for the long (solid) and short (dashed) graphene device, re- 
spectively. Three different values of the electron temperature 
in the electrodes were considered: 150 K, 300 K and 450 K. 
The dotted lines illustrate the I/ZcbT slope for the different 
temperatures. 

other hand, there is a strong temperature dependence, 
and in this case the electron transport is dominated 
by thermionic emission. The dotted lines illustrate the 
slope expected for thermionic emission. We see 
that in the gate voltage range from -0.25 to -0.75 V, the 
I-V characterics follow these slopes well. 

Fig. [6] also shows the electrostatic profile through the 
device. We see that the gate potential is almost per- 
fectly screened by the graphene ribbon, i.e. the gate po- 
tential does not penetrate through the central ribbon. 
This means that for a layered structure, only the first 
layer would be strongly affected by the back-gate. This 
has some implications also for gated nanotube devices. 
In such a device, only the atoms facing the gate elec- 
trode will be strongly influenced, and this explains why 
in Ref. 26 we found that the transport in the device was 
dominated by tunneling even though the nanotube was 
110 A long, and thus longer than the graphene junctions 
studied in this paper. Thus, to obtain efficient gating 
of a nanotube, the gate electrode must wrap around the 
tube. 



VI. CONCLUSIONS 

In this paper we have introduced a new semi-empirical 
model for electron transport in nano-devices. The model 
is based on the Extended Hiickel method that extends 
the work by Zahid et al}'^ to give a more complete de- 
scription of the electrostatic interactions in the device. In 
particular, the position of the electrode Fermi level and 
the charge transfer between the contacts and the device 
are calculated self-consistently. 

Compared to DFT-based transport methods, the main 
advantage of our new method is that it is computationally 
less expensive, as well as having the option of adjusting 
parameters to reproduce experimental data or computa- 
tionally very demanding many-body electronic structure 
methods. 

The model includes a self-consistent Hartree potential 
which takes into account the effect of an external bias as 
well as continuum dielectric regions and external electro- 
static gates. 

We used the model to study a Tour wire between gold 
electrodes, and found that the voltage drop in the device 
compares well with ah initio results, while the calculated 
current-voltage characteristics qualitatively agree better 
with experimental findings than the corresponding DFT- 
LDA results do. 

We also considered a graphene nano-transistor, and 
our study illustrated how the transport mechanism 
changes from tunnelling to thermionic emission as the 
device is made longer. 

These applications show that the new method can give 
an accurate description of a broad range of nano-scale 
devices. With its favorable computational speed, it is a 
good complement to ah initio-hdised transport methods. 
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